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Abstract 

We discuss the optimization of a neutrino factory for large sin^ 2613, where we assume minimum effort on the 
accelerator side. This implies that we use low muon energies for the price of an optimized detection system. We 
demonstrate that such a neutrino factory performs excellent if combined with the electron neutrino appearance 
channel. Instead of the platinum channel operated with the muon neutrinos from the muon decays, we propose to 
use the initial superbeam from the decaying pions and kaons, which might be utilized at little extra effort. Since 
we assume out-of-phase bunches arriving at the same detector, we do not require electron charge identification. In 
addition, we can choose the proton energy such that we obtain a synergistic spectrum peaking at lower energies. 
We find that both the superbeam and the neutrino factory beam should used at the identical baseline to reduce 
matter density uncertainties, possibly with the same detector. This effectively makes the configuration a single 
experiment, which we call "neutrino factory superbeam". We demonstrate that this experiment outperforms a 
low-energy neutrino factory or a wide band beam alone beyond a simple addition of statistics. 
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1. Introduction 

Future neutrino oscillation facilities are primar- 
ily targeted towards the discovery of a non-zero 
sin^ 2013, the neutrino mass hierarchy determina- 
tion, and the measurement of leptonic CP vio- 
lation. All of these measurements depend heav- 
ily on the magnitude of sin^ 2^13 . Two examples 
for experiments under way, which will test the 
magnitude of sin^ 26*13, are the T2K superbeam 
(SB) project [1] and the Double Chooz reactor 
experiment [2]. From all such neutrino oscilla- 
tion experiments, one may expect a sin^ 2^13 sen- 
sitivity limit of about 0.01 in a decade from now 
(see, e.g., Ref. [3]), i.e., sin^ 26*13 can either be 
constrained to be smaller than this value, or it 
will be discovered until then. Since the discovery 
of sin^ 26*13 affects the experimental strategy how 
to measure the mass hierarchy and CP violation, 
this "branching point" [4] will, at the latest, tell 
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US how to proceed. 

If sin^ 2013 turns out to be small, one may op- 
timize for a maximum sin^ 2613 reach of all mea- 
surements. It seems that a neutrino factory [5-8] 
(NF) or a higher gamma beta beam [9, 10] could 
be the most promising alternative; for the op- 
timization discussions, see, e.g., Refs. [11, 12]. 
For large sin^ 2^13, the measurement far away 
from the oscillation maximum and matter density 
uncertainties affect the NF performance [8, 13]. 
Therefore, other alternatives may outperform a 
NF, such as SB upgrades or beta beams (see, e.g., 
Refs. [11,14] for a direct comparison). 

Recently, potential improvements in the muon 
charge identification (CID) at low energies (see, 
e.g., Ref. [15]) have brought back a NF into the 
discussion for large sin^20i3. In Ref. [16], a 
muon energy of 4.12 GeV was proposed at base- 
lines of 1 000 to 1 500 km, an option similar to 
a wide band beam in terms of neutrino energies 
and baseline. The relatively low muon energy al- 
lows for a significant reduction of the effort on 
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the accelerator side, which means that this option 
might indeed be competitive to SB upgrades, etc.. 
Other alternatives to improve the NF for large 
sin^ 26'i3 are the use of the magic baseline [12, 17], 
or the platinum channel — > Ve (with the v^^s 
coming from the muon decays), which requires 
electron detection with CID [12]. In this way, 
the impact of matter density uncertainties can be 
significantly reduced as well. 

In this letter, we discuss the optimization of a 
low-energy NF. We compare such an optimized 
experiment to wide band beams and beta beams. 
Furthermore, we propose to use the neutrinos 
from the secondary pion/kaon beam in addition 
to the NF beam. This SB, possibly targeted to- 
wards the same detector, will provide similar in- 
formation to the platinum channel. However, it 
has two major advantages: First, the (at least 
in a magnetized iron calorimeter) very difficult 
electron CID is not required if operated out-of- 
phase with the NF beam. Second, the proton 
energy can be chosen such that the event rates 
are peaking at lower energies, where the NF has 
less events. We discuss the requirements for such 
a hybrid configuration, which we call "neutrino 
factory superbeam" (NF-SB), and we show the 
physics potential. 

2. Setup and Experimental Requirements 

The primary information from a NF comes 
from the "golden" Ve — > t^fj. appearance chan- 
nel [7]|1 The two major challenges for large 
sin^ 2^13 are the matter density uncertainties and 
the lack of events at low energies, where the inter- 
esting higher-order oscillation effects are present. 
As demonstrated in Ref. [12], complementary in- 
formation from the T-conjugated i'^ —^ V(, ap- 
pearance channel (platinum channel) would help, 
partly because of the correlated matter effect. 
However, CID for electrons in a magnetized iron 
calorimeter appears to be difficult because the 
electrons produce showers early (though it may 
not be impossible for low neutrino energies). In 
this letter, we propose to use the SB from the pion 
and kaon decays instead of the platinum channel. 

^Note we will also use antineutrinos in all of the discussed 
channels. 



The schematics for such a setup can be found in 
Figure [TJ We require that it be operated out- 
of-phase with the NF beam, which means that 
CID will not be required for this channel. In ad- 
dition, the proton energy can be adjusted within 
the range acceptable for the NF|f| such that a spec- 
trum peaking at lower energies can be obtained. 
As one can see from Figure [U there are two ma- 
jor requirements for such a setup: First, the ge- 
ometry has to be chosen such that the tertiary 
and quarternary neutrino beams point towards 
the same detector (not to scale). Second, the tar- 
get station has to be able to separate and focus 
the pion and kaon charges in a way that both ex- 
periments can be operated simultaneously. From 
Figure [U this requirement may be acceptable if 
we only use half of the luminosity of an exclusive 
NF (because only one polarity can be operated at 
the same time). However, one may speculate that 
a more sophisticated target station might recycle 
a fraction of the muons from the SB mode. Note 
that, for the NF plus SB, we have now three phe- 
nomenological degrees of freedom: Ep, E^, and 
L. 

The main prerequisite for a low-energy NF 
is an improved detection system with sufficient 
golden channel efficiencies at low energies. These 
efficiencies with a threshold possibly as low as 
0.5 GeV might be obtained by cuts leading to a 
higher charge mis-identification background, or 
a more refined detector than a magnetized iron 
calorimeter (such as a NOi^A-like detector or a 
liquid argon detector with a large magnetic vol- 
ume). For large sin^ 26*13, such a background will 
not be very important as long as the background 
level is < sin^ 26*13. We follow Refs. [12,16] by 
using a low-threshold (500 MeV) with improved 
energy resolution option (15%/-\/E), which we 
call Golden*. We us 50% efficiency for the golden 
channel requiring CID, and 90% efficiency for the 
disappearance channel without CID. Note that 
an even higher appearance efficiency may be ob- 
tained for the price of larger backgrounds. We 
use a total luminosity of 5 yr running time x 10^^ 
useful muon decays/year x 50 kt detector mass in 



*The acceptable range is relatively wide 8 GeV < Ep < 
30 GeV (see, e.g., Ref. [18]). ^ 
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Figure 1. Schematics of a NF-SB (not to scale). Our degrees of freedom are given by red/gray labels. 
Here we adopt the conservative point of view that only half the muons can be collected for the NF. 



each polarity (^"'", /i~ stored) if operated as NF 
alone, and half of that for the NF-SB. We assume 
that the number of useful muon decays/year ap- 
proximately requires a target power of 4 MW, but 
this number scales with non-trivial aspects of the 
muon cooling and acceleration system. For the 
comparison to the platinum channel, we again fol- 
low the optimistic choice in Ref. [12] using 40% 
efficiency, no upper threshold, and 1% charge mis- 
identification (c/., e.g., Ref. [19] for liquid argon). 

For the SB operation, we have tested a 
MiniBOONE-like wide band beam {Ep = 8GeV) 
and an AGS-like wide band beam {Ep — 28 GeV). 
Since the MiniBOONE-like beam has turned out 
too have too low neutrino energies for our ap- 
plication, we focus on the AGS-like beam in the 
following (see Refs. [14,20,21] for details on the 
simulation) . Note that a considerably higher pro- 
ton energy may lead to a significant loss of effi- 
ciency for the NF. We assume an electron detec- 
tion efficiency of 80% without CID (in Golden*) 
and that the experiments is subject to the same 
NC background as the muons. This assumption 
is probably too optimistic, but the intrinsic beam 
background is larger anyway. Our systematical 
errors are chosen to be 5% except for the signal 
normalization errors in the NF (2.5%). We com- 
pare in Figure [5] the event rate spectra from the 
three potentially used appearance channels (NF, 
SB, platinum). The peak of the platinum channel 
is at higher energies than the peak of the NF be- 
cause of the corresponding spectra from the muon 
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Figure 2. Appearance event rates as a function 
of the reconstructed neutrino energy. The fig- 
ures compare the NF (golden), SB, and platinum 
channel appearance. The curves are computed 
for E^ = b GeV, Ep = 28 GeV, L = 1 250 km, 
sin^ 26*13 = 0.1, (5cp = 0, WDot = 50 kt, 2.510^1 
useful muon decays, and Ptargot = 4 MW. 



decays (the electron neutrinos peak at lower en- 
ergies than the muon neutrinos). On the other 
hand, Ep is chosen such that the SB peaks at 
lower energies, where oscillation effects are more 
prominent (c/., right panel). Note that this figure 
assumes the muon decay numbers of the NF-SB, 
whereas an exclusive NF operation will allow for 
the double luminosity. 

For the simulation, we use the GLoBES soft- 
ware [22, 23]. For the oscillation parameters, 
we use sin^ 9i2 = 0.3, sin^ 023 = 0.5, Am^j = 
7.9 • 10-5 eV2, Ato2^ = 2.5 • IQ-^eV^, and a 
normal mass hierarchy unless stated otherwise 



4 




Figure 3. Fraction of Sep for which CP violation can be measured as a function of L and (sin^ 20i3 = 
0.1, 3tT confidence level). The left panel is for a NF alone with twice the luminosity as in the right 
panel. The right panel is for a NF plus SB combination with a fixed proton energy of Ep = 28 GeV. The 
contours are spaced by 0.02, where the thick curves correspond to a fraction of Sqp G {0.8,0.7,0.6, . . .}. 
The optima are marked by the diamonds. The luminosity "£" stands for 0.5 10^^ useful muon decays/yr. 



(see, e.g., Ref. [24]). In addition, we assume that 
Ato2i, 012, and the matter density will be known 
at the level of 5%. 

3. Optimization for Large sin^ 20i3 

As far as the optimizations of a low-energy NF 
and the NF-SB are concerned, let us first of all 
focus on the NF alone. In Figure [31 left panel, 
we show the fraction of (true) Sep for which CP 
violation can be measured as a function of L and 
Ef^ (for sin^ 2013 = 0.1). At the optimum (dia- 
mond), one can establish CP violation for 81% 
of all possible values of Sqp- One has close-to- 
optimum performance for L ^ 800 km — 1 300 km 
and Efj^ > 4 GeV, which is consistent with the 
configuration discussed in Ref. [16]. For the mass 
hierarchy, we find that L > 800 km to estab- 
lish the mass hierarchy for any Sqp (the E'^- 
dependence is moderate). In addition, note that 
for such a large sin^ 26*13, none of the discussed 
facilities has a problem to establish a non-zero 
sin^ 26*13. 

In Figure [3l right panel, we show the same op- 
timization for NF-SB with a fixed proton energy 
of Eri — 28 GeV and half the useful number of 



muon decays from the last example. At the op- 
timum (diamond), one can establish CP viola- 
tion for 85% of all possible values of Sep- The 
optimal baseline range depends on E^. For in- 
stance, one has close-to-optimum performance for 
L ~ 800 km - 1 500 km for E^, > 4 GeV. For the 
mass hierarchy, shorter baselines L > 500 km are, 
in principle, sufficient. This effect has been dis- 
cussed in Ref. [25]. Since the NF luminosity in 
this hybrid is only half of the one for the NF 
alone, we conclude that we have a identified a 
synergy between the NF and the SB which goes 
beyond the simple addition of statistics (for a def- 
inition of synergy, see Ref. [26]). Note that the 
SB alone operated at double luminosity approxi- 
mately compares to the low-energy NF alone (for 
sin^ 2^13 = 0.1 and the same detector). 

Now one may argue that it may be better to 
use a combination with two different baselines in 
the fashion of Ref. [27] (NF plus SB). Therefore, 
we show in Figure 3] the two-baseline optimiza- 
tion of the NF and SB with uncorrelated matter 
effect and a fixed muon energy of E^ = 5 GeV. 
The optimum turns out to be close to the same- 
baseline diagonal. However, it is somewhat worse 
than for the correlated matter density case above. 
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Figure 4. Two-baseline optimization of the NF 
plus SB with uncorrelated matter effect: Fraction 
of Sep for which CP violation can be discovered 
as a function of Lnf and Lsb (3ct). The muon 
energy is = 5 GeV, and sin^ 2^13 = 0.1. 



Since a option requiring only one site and, possi- 
bly, only one detector is much more cost efficient, 
we will not discuss the two-baseline case any- 
more. We have, in addition, tested a 500 kt wa- 
ter Cherenkov detector for the SB part together 
with the 50 kt Golden*. More or less, we obtain 
the same qualitative results with the difference 
of a much higher fraction of 6cp- CP violation 
can be measured for about 91% of all values of 
Sep (for sin^ 2013 = 0.1). Of course, this option 
would require an additional very large detector. 
Furthermore, we have checked that our optimiza- 
tion results to not change for somewhat smaller 
sin^ 2^13 = 0.05. In the following section, we will 
use these standard setups from the optimization 
discussion above as given in Table [T] 

4. Comparison of Facilities 

We compare in Figure [5] the different se- 
tups from Table [T] for different values of (large) 
sin^ 2^13. Obviously, the NF-SB outperforms any 
of the chosen facilities for large enough sin^ 26'i3 
except for the beta beam. In particular, it is sig- 
nificantly better than low-E NF (with or without 
platinum channel), and slightly better than High- 
E NF2B (except for sin^ 26*13 = 0.01). Note that 



these setups are operated at the double muon de- 
cay rate compared to the NF-SB, which means 
that the information from NF and SB is very com- 
plementary in terms of channel and spectrum. As 
far as the beta beam is concerned, we find that 
it has the best absolute performance. However, 
given the relatively high gamma and the challeng- 
ing ion decay rates, it is difficult to judge how 
this setup compares in terms of feasibility and ef- 
fort. The comparison to the wide band beam (in 
this case with a 500 kt water Cherenkov detector 
compared to the 50 kt Golden* used for the NF- 
SB hybrid) depends very much on sin^ 2^13: For 
very large sin^ 26*13 — 0.1, the wide band beam 
might be the natural choice, because it is based 
on known technology. However, for a somewhat 
smaller sin^20i3, the NF-SB is significantly bet- 
ter. 

We have also studied further potential improve- 
ments of the NF-SB. Very little physics reach can 
be gained by the additional use of the magic base- 
line or platinum channel, or by better knowl- 
edge on the matter density profile, since the un- 
derlying correlations are already resolved. The 
most significant improvement can be obtained by 
luminosity upgrades, such as by using a large 
water Cherenkov detector (for the SB) in addi- 
tion to Golden* at the same site. In this case, 
even the beta beam can be outperformed for 
very large sin^ 2^13. Note that one could imagine 
different staged scenarios with this observation, 
such as WBBGoiden* followed by NF-SB. We ex- 
pect that further significant improvements could 
be obtained by better absolute efficiencies for the 
muon neutrino appearance channel (for the price 
of larger backgrounds), by recycling a fraction of 
the muons produced in the SB channel for the NF, 
or by a further optimization of proton energy and 
horn characterizing the SB. 

5. Summary and Discussion 

We have studied a low-energy neutrino factory 
(NF) in combination with a superbeam (SB) pro- 
duced in the decay chain of the same protons, 
which we have called "neutrino factory super- 
beam". We have demonstrated that this config- 
uration outperforms a low-energy NF alone be- 
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Setup name 


Description/Parameters 


Dot. 




Ref. 


WBBwc 


Ep = 28GeV, Pt^g^t =4MW 


WC 


500 kt 


[20] 


Low-E NF 


E^, = 5 GeV, L = 900 km, lO^^ umd 


Golden* 


50 kt 


this 


Low-E Nf&iPt 


Same as Low-E NF, plus platinum channel 


Golden* 


50 kt 


this 


High-E NF2B 


Ef, = 20 GeV, L = 4 000 km + 7 500 km, lO^i umd 


Golden* 


2 X 50 kt 


[12] 


NF-SB 


Ef,^5 GeV, L = 1 250 km, 0.5 lO^i umd. 


Golden* 


50 kt 


this 




Ep - 28 GeV, Ptargct = 4MW (SB) 








BB350 


7 = 350, L = 730 km, 5.8 • lO^^ useful ^He 


WC 


500 kt 


[10] 



decays/yr, 2.2 • 10^^ useful ^*Ne decays/yr 
Table 1 

Setups used for the comparison. The running time is five years for each polarity for all setups. The abbre- 
viation "umd" stands for "useful muon decays per year" , the abbreviation "WC" for a water Cherenkov 
detector, and the label Golden* for our optimized detector (see description in Section [5] and Ref. [12]). 
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Figure 5. Comparison of the fraction of (true) (5cp for which CP violation can be discovered (3cr) for the 
different facilities listed in Table [TJ The different panels correspond to different values of sin^20i3, as 
given in the captions. 



yond a simple addition of statistics even if one 
can use the platinum channel, as well as it out- 
performs a wide band beam with the same proton 
energy. The reason for this excellent performance 
is the SB spectrum peaking at lower energies, and 
the information from the T-conjugated —^ 
channel at the same baseline (i.e., with correlated 
matter effect). 

As far as the technical requirements for such 
a configuration are concerned, the same detector 
might be used at the same baseline. As for any 
low-energy neutrino factory, sufficiently good ef- 
ficiencies at low energies are required. The extra 



effort mainly reduces to a more sophisticated tar- 
get station, which has to separate the pion and 
kaon charges and utilize them for the NF and 
SB. In addition, we have assumed that the NF 
and SB bunches arrive out-of-phase at the (pos- 
sibly same) detector. For the electron neutrino 
detection, we have required a relatively high ef- 
ficiency without CID, which is different from the 
platinum channel (where the requirement of CID 
significantly reduces the efficiency). As far as sys- 
tematics is concerned, we expect this experiment 
to be of low risk: The NF systematics has been 
assumed very conservatively, and the wide band 
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beam is rather robust with respect to systemat- 
ics [14]. This might be an important decision cri- 
terion in comparison to narrow band beams. 

From the optimization, we have found L ~ 
800 km - 1 500 km and E^>A GeV to be close- 
to-optional, which, for instance, corresponds to 
Fermilab to Homestake or Henderson. For the 
proton energy, we have tested two configurations 
at the lower end upper end of the admissible range 
for a NF, and we have chosen Ep = 28 GeV. 
This means that we obtain a completely differ- 
ent configuration from a high-energy NF opti- 
mized for small sin^ 2^13 [12] in terms of base- 
fines, muon energy, and target station. For the 
detector, the need for high efficiencies at low ener- 
gies is independent of the NF configuration (but 
not a prerequisite for a high-energy NF). How- 
ever, higher backgrounds may be admissible for 
the low-energy NF, which means that the over- 
all efficiency could be higher. Although a high- 
energy NF with two baselines (and possibly plat- 
inum channel) might be used for large sin^ 2^13 
as well, the effort might be considerably higher. 
The same, in principle, applies to a higher gamma 
beta beam. We conclude that a reasonable deci- 
sion can be made when the branching point sen- 
sitivity sin^ 2^13 ~ 0.01 is reached. Before that 
point, one may want to develop two different NF 
approaches: One optimized for sin^ 2^13 reach op- 
erated as a discovery machine, the other opti- 
mized for large sin^ 2^13 operated as a precision 
instrument. 
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